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ABSTRACT 


The knowledge of the ocean environment, especially its 
littoral/coastal waters, is important for current and 
future naval operations. In particular, an accurate 
description of the mesoscale variability of the surface and 
subsurface currents is essential for optimizing Navy 
missions. 

This study determines that different types of 
topographic smoothing and the addition of a volume 
constraint to a sigma-coordinate model significantly 
influences the generation, evolution, and maintenance of 
not only the surface and subsurface currents, but also 
other features (upwelling, meanders, eddies, filaments, 
Mediterranean Outflow, and Meddles) in the Northern Canary 
Current System (NCCS) . The NCCS is chosen for this study 
for its classical upwelling and equatorward surface current 
as well as the unique impact of the Mediterranean Outflow. 

The sound speed structure of these features is 
examined to establish which regions of the NCCS experience 
the largest changes in sound speed and the most intense 
gradients. This study demonstrates that features of classic 
littoral/coastal eastern boundary current systems, like the 
NCCS, have an important impact on the sound speed structure 
which can significantly impact Navy sonar operations. 
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I. 


INTRODUCTION 


An accurate description of the mesoscale variability 
of surface and subsurface currents is essential for 
optimizing Navy missions. Understanding the mesoscale 
dynamics of the world's oceans is crucial to building data 
bases and developing operational forecasting models of the 
sea environment for defense-related operations. 

As naval operations continue to occur in coastal and 
littoral/coastal regions, accurate predictions of shallow 
water acoustics are essential. Although the sound speed 
structure in the open ocean is well known, the complex 
dynamics of the littoral/coastal regions can complicate 
acoustic prediction near the coast. Identifying the 
influence of dynamical processes such as upwelling, eddies, 
filaments, and meanders on sound speed structure can 
provide better ocean models for acoustic predictions in 
littoral/coastal regions. 

This study determines if different types of 
topographic smoothing (e.g., Gaussian and direct iterative 
methods) and the addition of a volume constraint to a 
sigma-coordinate model (in this case, the Princeton Ocean 
Model) can significantly influence the generation, 
evolution, and maintenance of currents, upwelling, 
meanders, eddies, filaments, Mediterranean Outflow, and 
Meddles in the Northern Canary Current System (NCCS). In 
this section, the descriptive, dynamic and modeling 
background of the NCCS region is reviewed. 
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A. NORTHERN CANARY CURRENT SYSTEM (NCCS) 

The NCCS forms the eastern limb of the subtropical 
gyre of the North Atlantic. Flowing equatorward off the 
west coasts of the Iberian Peninsula (IP) and Morocco, the 
NCCS contains classical features of an eastern boundary 
current (EEC) system, as well as regional anomalies. These 
features are wind-driven and are characterized by a 
climatologically weak (<10 cm/s), broad (-1000 km wide), 
cold surface flow towards the equator. The surface current 
flows in the direction of the prevailing winds with a 
poleward undercurrent below. The EEC thermocline is 
typically shallow (<30 m depth) and the region has high 
biological production due to vast regional upwelling 
(Parrish et al., 1983). While the NCCS contains a number of 

these classical features, the region also includes notable 
anomalies. 

The Canary Current (CC) extends south along the 
western coasts of both the IP and northwest Africa from 
~10°N to ~45°N. Like other EECs, the NCCS is a relatively 
weak (-10-30 cm/s), broad (-1000 km), shallow (-500 m) 
surface current that flows equatorward year-round (Wooster 
et al., 1976; Eatteen et al. , 2000). Its surface waters are 

relatively cool due to upwelled water from the coast 
(Mittelstaedt, 1991). 

Characteristic of EECs, a subsurface poleward 
undercurrent flows below the surface CC over the 
continental slope and is relatively narrow (-10-40 km) and 
weak (-2-10 cm/s) . Varying with season and latitude, the 
undercurrent is strongest between depths of -100 m and -600 
m. In winter, the CC shoals to the north near Cabo da Roca 
and forms a second surface current commonly referred to as 
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the Iberian Current (IC) (Haynes and Barton, 1990) . The IC 
is a relatively narrow (~25-40 km), weak (~20-30 cm/s), 
seasonal current found trapped near the coast against the 
shelf break (Fiuza, 1980; Frouin et al. , 1990; Haynes and 

Barton, 1990) and is usually found north of Cabo da Roca 

but can be seen as far south as Cabo de Sao Vicente 

(Batteen et al., 2000). 

The NCCS is unmatched amongst the subtropical eastern 
boundary currents of the world in that it contains the 

unique influence of Mediterranean Outflow (MO). After 

leaving the Strait of Gibraltar, the MO descends close to 
the continental slope in the Gulf of Cadiz as a narrow 
gravity current. This salty MO plume is deflected to the 
right due to the Coriolis acceleration (Bower et al., 1997) 

and dilutes, thickens, and becomes vertically 
differentiated into two distinct cores as it follows the 
northern rim of the Gulf of Cadiz (lorga and Lozier, 1999; 
Bower et al. , 1997) . At ~7°W in the Gulf of Cadiz, the two 

cores exist at depths of ~600-900 m and ~1100-1200 m. Both 
cores continue to flow westward along the southern coast of 
Spain and turn poleward around Cabo de Sao Vicente (Amber 
and Howe, 1979; lorga and Lozier, 1999) . A third, 
shallower, poleward core of MO has also been traced from 
the Strait of Gibraltar northward to ~38.5°N off western 
Portugal (Amber, 1983). Additionally, climatological 
cyclonic circulation in the southwestern Gulf of Cadiz 
spreads the salty MO south to ~34°N (lorga and Lozier, 
1999) . 

Typical of other EBCs, the NCCS contains highly 
energetic mesoscale features, including jet-like surface 
currents, meanders, eddies and filaments (Johnson and 
Stevens, 2000). Satellite observations during periods of 
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upwelling-favorable winds have shown nearshore upwelling 
with several narrow filaments of cooler water extending off 
the coast of the IP (Fiuza and Sousa, 1989) and Cape Ghir 
in northwest Africa (Van Camp et al. , 1991; Hagen et al. , 
1996). Upwelling filaments often extend ~80-150 km offshore 
with alongshore spacing of ~80-100 km between filaments 
(Haynes et al. , 1993) and terminate with dipole eddy pairs 
(Fiuza et al., 1982; Barton, 1998). 

Anticyclonic and cyclonic pairs of mesoscale eddies 
have been observed off the coast of the IP on the order of 
100 km (Fiuza, 1984; Stammer et al., 1991). These mesoscale 
features have also been observed during periods of 
upwelling-favorable winds and are located near coastline 
irregularities, such as capes. These observations provide 
evidence that wind forcing, along with coastline 
irregularities, appears to be an important mechanism in the 
formation and continuance of many of the mesoscale features 
found in the NCCS (Batteen et al., 2000). 

Unique to the NCCS is the generation of subsurface, 
anticyclonic, submesoscale eddies, called Meddles, which 
contribute to the maintenance of the Mid-Atlantic Salt 
Tongue (e.g.. Arm! et al. , 1989) . Meddles form from the 
complex flow of waters entering and exiting through the 
Straits of Gibraltar; cool North Atlantic water flows into 
the Mediterranean, while at the same time warmer, denser 
and more saline water exits beneath the incoming flow. The 
high-salinity water MO flowing out of the Mediterranean Sea 
descends to mid depths in the density-stratified ocean, 
continues as a narrow jet along the Iberian continental 
slope, and intermittently detaches large-scale eddies, 
called Meddles. This is an important process because it 
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maintains the relatively high mean salinity of a major 
water mass in the North Atlantic Ocean. 

The primary generation region of Meddles is off Cabo 
da Roca and Cabo de Sao Vicente off southwest Portugal. 
Several different trajectories of Meddles have been 
observed, including a southwestward movement into the 
Canary Basin and westward translations south of the Azores 


(Richardson 

and Tychensky, 1998). 

It 

has been 

shown 

that 

bathymetry 

impacts the evolution 

of 

Meddles 

(Price 

and 


Baringer, 1994). The presence of sloping bathymetry along 
the continental margin, combined with the general 
circulation, generates a complex flow pattern along the 
slope regions. The degree to which the MO descends along 
the continental slope is sensitive to the steepness of the 
bottom slope it travels along. As the MO transits along the 
slope. Meddles can also spin off to form isolated eddies at 
abrupt changes in topography, e.g., as observed off Cabo de 
Sao Vicente (Serra et al., 2002; Bower et al., 1997). 

Anticyclonic vortexes associated with Meddles form 
through frictional flow along the IP. Transiting along the 
coast of Portugal at ~2-5 cm/s. Meddles propagate to the 
southwest due to anticyclonic shear and interaction with 
the southward flow of the stronger CC. Observations reveal 
that Meddles can contain a core of Mediterranean waters 
-200 to -1000 m thick centered at about -1000 m and are 
less stratified than the surrounding water. They are 
typically -20-100 km in diameter and have temperature and 
salinity anomalies of approximately -2°C and -1 psu. This 
density structure supports an energetic anticyclonic 
rotation with peak azimuthal velocities of about -30 cm/s 
and rotation periods of -2-7 days (Arm! et al., 1989; 

Hebert et al., 1990). 
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B. PREVIOUS MODELING STUDIES 

Previous modeling studies (e.g., Batteen et al., 2000) 
have focused on the driving mechanisms of complex mesoscale 
activities in the NCCS region, including upwelling, 
filaments, highly energetic eddies, and meanders. Martinho 
(2003) used the results of four numerical experiments of 
increasing complexity to investigate the roles of wind 
forcing, bottom topography and thermohaline gradients on 
the features of the NCCS. The results of this study showed 
that while wind forcing is the primary mechanism for 
generating classical EBC features, bottom topography and 
thermohaline gradients can play important roles in the 
generation, evolution, and maintenance of classical as well 
as unique features in the NCCS. 

In a more recent study, Guthrie (2006) investigated 
the effect of currents, wind and eddies on the sound speed 
structure of two different EBC models, the NCCS and the 
Leeuwin Current System. She found that the dynamics of EBCs 
can lead to large changes in sound speed (5-15 m/s) and 
distort the vertical sound speed profile. 

C. ORGANIZATION 


The objective of this study builds on these previous 
studies. A new volume constraint and different types of 
bottom topography to the model are input and investigated 
to determine their effect on generation, evolution, and 
maintenance of currents, upwelling, meanders, eddies, 
filaments, MO, and Meddles, as well as on the sound speed 
structure in the NCCS. Knowledge gained from this study 
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will allow Navy oceanographers to better understand the 
complex ocean dynamics observed in EBCs and more accurately 
predict sound speed variability in these regions. In turn, 
the Navy will be better equipped to support mine and anti¬ 
submarine warfare operations in these areas. 

The organization of this thesis is as follows: Chapter 
II describes the Princeton Ocean Model (POM) used 
throughout the study, including its data sets, 
preprocessing, and boundary conditions. In Chapter III, new 
numerical methods used to improve modeling capabilities of 
sigma coordinate bottom-following models are examined and 
compared. 

The results of several process-oriented experiments 
(e.g., see Table 1) are discussed in Chapter IV. These 
experiments, arranged in order of increasing complexity, 
are as follows: Experiment 1 isolates the effect of wind¬ 
forcing. Experiment 2 examines the additional role of 
topography by incorporating a Gaussian type of bottom 
topography. Experiment 3 investigates the additional role 
of thermohaline forcing. Einally, Experiment 4 examines the 
sensitivity of the type of bottom topography used by 
incorporating an iterative instead of Gaussian type of 
bottom topography. These experiments seek to isolate 
phenomena in the current system that is generated by the 
various types of bottom topography. Lastly, the effects of 
these dynamic features are explored on the sound speed 
structure. Eindings, conclusions and recommendations are 
presented in Chapter V. 
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Experiment 

Nvimber 

Annual Climatology 

Bottom 

Topography 

1 

Horizontally 

Averaged 

Flat 

2 

Horizontally 

Averaged 

Smoothed with a 

2-D Gaussian 

Filter Method 

3 

Full 

Smoothed with a 

2-D Gaussian 

Filter Method 

4 

Full 

Smoothed with 
the 1-D Direct 

Iterative 

Method 


Table 1. 


Summary of experimental design. 




II. MODEL DESCRIPTION 


The sigma-coordinate Princeton Ocean Model (POM), 
established by Blumberg and Mellor (1983, 1987), is well 
recognized and widely used in numerical oceanography. The 
POM is a three-dimensional (3-D), sigma coordinate 
(terrain-following), primitive equation ocean model with a 
free surface. One unique feature of this model is its 
embedded second order turbulence scheme, which yields 
realistic Ekman surface and bottom boundary layer flows and 
effectively suppresses two-grid point noise in the model 
output (Mellor and Yamada, 1982). 

Through the use of bottom-following sigma levels, the 
POM can realistically simulate processes associated with 
continental shelf flows and bottom boundary layer dynamics 
in local domains in coastal regions. Regions of complex 
topography, such as the NCOS, are simulated well with sigma 
coordinate models. POM has been used successfully to model 
coastal regions, where variations in bottom topography are 
important, and more recently, to model entire ocean basins 
(Ezer and Mellor, 1994, 1997) . 

Here the POM uses 21 vertical sigma levels which are 
most closely spaced near the surface and bottom. The 
vertical resolution is higher near the surface and the 
bottom in order to resolve both the surface boundary layer 
and the bottom boundary layer, which are important in 
coastal regions. To eliminate the time constraints for the 
vertical grid related to higher resolution near the 
surface, bottom and in shallow waters, an implicit vertical 
time differencing scheme is used. 
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The prognostic variables of the model are: potential 
temperature, salinity, density, three velocity components, 
surface elevation and turbulent kinetic energy (Blumberg 
and Mellor, 1987). The model has a split time step interval 
for external and internal modes. The external mode solves 
the equations for the vertically integrated momentum 
equations and also provides the sea surface and barotropic 
velocity components with a time step of six seconds. The 
internal mode solves the 3-D equations and has a time step 
of 180 seconds in the NCOS model. 

The horizontal diffusion is a Smagorinsky formulation 
(Smagorinsky et al., 1965) in which the horizontal 
viscosity coefficients depend on the grid size, the 
velocity gradients and a slope parameter coefficient. 
Consistent with other POM studies, a value of 0.2 is 
assigned to this coefficient (Ezer and Mellor, 1997) . This 
reduces unwanted diapycnal mixing over steep topography. 

All choices for processes and values used in the model 
are the default POM parameters unless specifically 
described otherwise. This formulation is documented in 
detail by Mellor (2004). 

A. DATA SETS 

The model domain (Figure la) encompasses the west 
coasts of Iberia and Morocco, extending from 30°N to 41.5°N 
and from 5°W to 16.5°W. The coastal boundaries of the model 
domain are closed, including the Strait of Gibraltar. The 
northern, western, and southern borders are open 
boundaries. 

Topographic data was collected from the Institute of 
Geophysics and Planetary Physics, University of California 
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San Diego (Sandwell and Smith, 1996) from a compilation of 
30 years of data soundings by ships. In vast regions, where 
there is sparse or no ship data, gravity estimates, derived 
from altimetry data, are used to interpolate soundings. The 
bottom topography has a horizontal of 2 minutes (1/30 of a 
degree). 

The model is initialized using annual climatological 
temperature and salinity values obtained from Levitus and 
Boyer (1994) and Levitus et al. (1994). This data set has a 
one-by-one degree horizontal resolution. Climatological 
wind fields are obtained from the European Centre for 
Medium Range Weather Forecasts (ECMWF) near-surface wind 
analyses (Trenberth et al., 1990). The wind data is 
formulated on a 2.5-by-2.5 degree grid. All experiments are 
performed on a beta plane with a realistic coastline. 

B. PRE-PROCESSING 

Pressure gradient force errors (PGFE) in a sigma 
coordinate model are inevitable due to sloping sigma 
levels. The main concern when using a terrain bottom¬ 
following ocean model is to reduce the PGFE. The highest 
resolution in the model domain is where the values of the 
slope parameter, defined as the absolute value of the ratio 
of the difference between two adjacent grid cell depths and 
their mean depth, are largest in both the latitude and 
longitude directions (Mellor et al. , 1998) . According to 
Mellor et al. (1998), the maximum suggested slope parameter 
to be used in POM sigma coordinate models is 0.2, since 
greater values can induce higher PGFEs. Minimizing the 
slope parameter is accomplished by constructing the grid 
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such that the closest grid spacing is located over the 
steepest regions of the continental slope. 

Over much of the model domain, the steepest bottom 
gradient still had a slope parameter greater than 0.2. 
Therefore, the raw topography had to be smoothed in order 
to minimize the PGFE. Two very different linear 
interpolation filters used to reduce the slope parameter 
are discussed in Chapter III. 

The annual temperature and salinity values were 
obtained from the World Ocean Atlas 1994 and interpolated 
from the one-by-one horizontal grid at standard depths to 
the varying horizontal resolution and 21 vertical sigma 
levels of the POM using a 3-D linear interpolation scheme. 
Due to the changes in the vertical levels for each 
different type of topography used, this was done separately 
for the three different topographies used in this study: 
flat, smoothed, and iterative topographies. 

Daily seasonal winds were averaged over time in order 
to form the annual, non-weighted, average wind vector 
field. The wind vectors were also interpolated for the 
horizontal resolution of the POM with a two-dimensional (2- 
D) linear interpolation scheme. Since wind stress forces 
the model, the annual wind vectors were then used to 
compute the wind stress. In all experiments, annual wind 
forcing is used. 

C. INITIALIZATION, FORCING AND BOUNDARY CONDITIONS 

As discussed above, POM was initialized with annual 
temperature and salinity values obtained from Levitus and 
Boyer (1994) and Levitus et al. (1994) . Since the model 
runs reached a quasi-equilibrium state in a relatively 


12 



short time (<26 days), zero salinity and temperature fluxes 
have been prescribed at the ocean surface. Experiments 1 
and 2 are initialized with horizontally averaged 
climatology which uses the same initial temperature and 
salinity over the entire model domain. Initial full 
temperature (Figure 2) and salinity (Figure 3) climatology 
of the NCOS region is used in Experiments 3 and 4. 
Temperature and salinity decrease roughly linearly from 
north to south by 3°C and 0.8 psu, respectively. 

The model was forced from rest with the annual ECMWF 
wind fields, which were interpolated for the model grid. 
This wind field (Figure 4) is northerly, towards the 
equator, along the entire domain and increases in magnitude 
towards the south. 

Currently, no general criteria are available that 
state which boundary conditions are best used for a 
specific model. Specifying the correct open boundary 
conditions (BCs) is vital for achieving realistic results. 
These realistic results have no reflections, clamping, 
spurious currents or numerical alteration of the total 
volume of water in the model. For free-surface models, such 
as the POM, one important criterion is that the BCs should 
be transparent to the waves. In this model, a gradient BC 
(Chapman, 1985), which allows geostrophic flow normal to 
the boundary, is used for the elevation. For inflow 
situations, the model is forced with annual temperature and 
salinity values obtained from Levitus and Boyer (1994) and 
Levitus et al. (1994), while for outflow situations, an 

advection scheme is used. 

For the baroclinic velocity components normal to the 
boundary, an explicit wave radiation scheme based on the 
Sommerfield radiation condition is used. For the barotropic 
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velocity components, a radiation condition, the Flatter 
radiation plus Roed local solution was used. Palma and 
Matano (2000) revised the Flatter (1976) radiation 
condition using the Roed and Smedstad (1984) technique, 
which breaks the velocity into both global and local 
components. This radiation condition allows wave energy to 
radiate across the boundary while conserving mass inside 
the model domain. 

Even though the topography field produces a closed 
Strait of Gibraltar, the MO is represented as an initial 
condition by the annual temperature and salinity values. 
Since the model runs a quasi-equilibrium state, in a 
relatively short manner (~26 days, not shown), the ongoing 
addition of heat and salinity is not required to maintain 
the MO. 

D. SOUND SPEED CALCULATION 

Sound speed is a complex function of temperature, 
salinity, and pressure (Fofonoff and Millard, 1983) . In 
this study, the UNESCO international algorithm (equation) 
is used to calculate the speed of sound in seawater. This 
sound speed equation is based on the equation of Chen and 
Millero (1977) with coefficients recalculated by Wong and 
Zhu (1995). This equation is far more complicated than most 
others in use as pressure is used as a variable while 
simpler forms use depth. 

Note that a number of other equations have been 
published for calculating the sound speed in seawater based 
on laboratory measurements (e.g., Wilson, 1960; Del Grosso, 
1974; Coppens, 1981; and Mackenzie, 1981). Each of these is 
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different and there is much debate over the accuracy and 
range of applicability of each. In particular, Dushaw et 
al. (1993) has shown that there is a difference of 0.6 m/s 
at 5000 dbar between these equations (e.g., the Del Grosso 
equation and the Chen and Millero equation) . Since the 
model results in these experiments extend to 2500 m depth, 
the Chen and Millero equation was chosen due to its 
international acceptance and greater range of applicability 
(i.e., applicable to shallow and deep depth as opposed to 
only a limited range for simpler equations). 
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Ill. MODEL IMPROVEMENTS 


A. BOTTOM TOPOGRAPHY INCORPORATION 

Velocity errors induced by the pressure gradient force 
are unavoidable in 3-D sigma-coordinate models, such as the 
POM. PGFEs are produced when computing the horizontal 
pressure gradient near steep topography in sigma-coordinate 
oceanic models (Mellor et al. , 1994, 1998). Significant 
error is caused by the separation of the horizontal 
pressure gradient term into two parts and the subsequent 
cancellation of the two terms in the presence of steeply 
sloping topography, which results in large truncation 
errors in the transformation from z-coordinates to sigma- 
coordinates (Luo et al., 2002). 

The main concern when using a terrain bottom-following 
ocean model is to reduce the pressure gradient force error 
(PGFE). One approach used to reduce the error is to 
decrease the slope parameter. Here, two interpolation 
filter techniques are discussed. One method to reduce the 
slope parameter is with the traditional 2-D smoothing with 
Gaussian filter, provided by the POM. The other method is a 
one-dimensional (1-D) robust direct iterative technique. 
This method has been shown to have the unique advantage of 
maintaining coastline irregularities, continental shelves, 
and relative maxima such as seamounts and islands (Martinho 
and Batteen, 2006) . 

As listed in Table 1, Experiment 1 will use flat 
bottom topography with a constant depth of 2500 m. For both 
Experiments 2 and 3, smoothed bottom topography will be 
incorporated. The raw topography (Figure lb) is smoothed 
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with a linear, 2-D low-pass Gaussian filter in order to 
decrease the pressure gradient error to the limit of 0.2 at 
all grid points. The new depth of each point calculated 
with this filter has a non-weighted average of 15-by-15 
points surrounding the point. Depths greater than 2500 m 
are reassigned to the maximum depth of 2500 m. Land is 
assigned a depth of 20 m to avoid division by zero in the 
model. The Strait of Gibraltar is closed. Figure Ic shows 
the topography after the Gaussian smoothing which is used 
in Experiments 2 and 3. 

For Experiment 4, raw topography is smoothed with an 
iterative, 1-D technique (Martinho and Batteen, 2006) which 
also satisfies the slope criterion. This is an iterative 
process since a change in the topography necessary to 
reduce the slope parameter in one direction may alter the 
slope parameter in the perpendicular direction to values 
greater than 0.2. First, the 2-D Gaussian filter is applied 
to the coastline to avoid altering the adjusted slope 
parameters. Then, the slope parameter is calculated along 
each grid line in one, initial direction over the domain. 
Any two adjacent grid points with a slope parameter greater 
than 0.2 is adjusted and the line recalculated. After all 
grid lines in the initial direction are adjusted, the 
topography matrix is rotated 90 degrees. Each process is 
repeated until the topography is adjusted for all 
directions (rotated 360°). After the smoothing, land is 
assigned a depth of 20 m and depths greater than 2500 m are 
reassigned to 2500 m in order to avoid division by zero in 
the model. Again, the Strait of Gibraltar is closed. Figure 
Id shows the topography after the iterative smoothing, 
which is used in Experiment 4. 
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For comparison, 3-D plots of the raw, Gaussian, and 
iterative topographies for the NCOS are shown in Figures 
5a, 5b and 5c, respectively. There is a recognizable 
difference between the Gaussian (Figure 5b) and raw 
topography (Figure 5a). A significant widening of the 
topographic features can be seen in the Gaussian smoothed 
topography. There is a clear deepening, or "sinking", of 
seamounts and islands and a significant smoothing of the 
coastline irregularities. For the continental slope and 
rise regions, there is a relatively large change in depth. 

In contrast, no noticeable differences are seen 
between the raw (Figure 5a) and the iterative (Figure 5c) 
topography. The minimum depth of the seamounts is 
maintained which better preserves the continental shelf and 
rise and coastline irregularities. The iterative method 
only changes the topography in areas where there are 
relatively high slopes and shallow depths: near the upper 
continental slope and around seamounts. This slight 
widening of the most prominent seamounts is shown in Figure 
5d. 

B. INCORPORATION OF A VOLUME CONSTRAINT 

An advanced volume constraint (VC) based on research 
performed by Marchesiello et al. (2001) is introduced to the 
POM as a subroutine. The VC, an adaptive algorithm to the 
model, is based on global correction of normal barotropic 
velocities. 

Free-surface models, such as the POM, can lose or gain 
water through the boundaries. According to Marchesiello et 
al., (2001), it is necessary to apply an artificial volume 
constraint to the model since the algorithm for the 
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barotropic velocities is non-volume conservative. In order 
to conserve volume, the total volume of water transported 
through the boundaries is calculated at every external time 
step. If there is a net volume outward or inward, a 
correction is applied at every boundary point to obtain a 
net change in volume equal to zero. Since the correction is 
always small and is applied at every open boundary 
condition point, there should be no marked change in the 
flow structures (Marchesiello et al., 2001). 

To test whether the VC is an improvement to the NCCS 
model, two runs were made: one with a VC and one without a 
VC. Even by day 1 (Figures 6a and 7a), discernible 
differences are noted. For example, the surface flow is 
stronger to the south in the run without the VC (Figure 7a) 
compared to the run with the VC (Figure 6a) . In time, the 
model without the VC (Figure 7b) continues to diverge from 
the model with the VC (Figure 6b). In particular, the model 
run with the VC (Figure 6b) adjusts faster and shows more 
realistic results, i.e., equatorward flow to the south 
(consistent with the wind forcing) with no spurious 

poleward inflow from the southern boundary of the model. 

Consistent with the findings of Marchesiello et al. 
(2001), adding the VC improved model stability, helped 
avoid accumulating errors on the long term, and improved 

the overall performance of the open BCs in the POM. 

Compared to previous studies on the NCCS which had no VC 

(i.e., Martinho, 2001; Guthrie, 2006), this improvement of 
adding the VC drastically reduced the loss of the volume at 
the boundaries to insignificant levels throughout the model 
domain. As a result, the incorporation of the VC is used in 
all the experiments listed in Table 1 for the NCCS. 
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IV. RESULTS 


A. EXPERIMENT ONE 

In Experiment 1 (see Table 1), the model was 
initialized with the horizontally averaged annual 
climatological temperatures and salinities, which removes 
all thermohaline forcing. A flat bottom with a constant 
depth of 2500 m eliminated any topographic effects. A 
realistic coastline was used and the model was forced from 
annually climatological winds. Using these initial 
conditions, the model was run for 90 days. 

1. Surface Currents and Upwelling 

By day 30 (Figure 8a), as expected, there is an 
equatorward surface current along the entire coastal 
region, with stronger surface currents found in the region 
between Cape Bedouzza and Cape Ghir (see Figure la for 
geographical locations) . The stronger winds at the 
southward end of the model cause cooler temperatures 
associated with relatively strong upwelling along the 
entire coastline. By day 45 (Figure 8b), the offshore 
extent of the coastal upwelled waters has increased, 
particularly off the coast of Cape Ghir, where the wind 
stress is initially stronger. By day 75 (Figure 8d) , the 
coldest waters are found off Cape Ghir and Cabo da Roca, 
and the upwelling extent continues to broaden along the 
entire coastline. This continuing westward propagation of 
upwelled water with time is consistent with the results of 
McCreary et al. (1987) and Batteen et al. (1989) . They 

showed that, due to the beta effect, the surface coastal 
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jet does not necessarily have to be confined to within a 
Rossby radius of deformation of the coast. 

2. Undercurrent 

Below the surface, a poleward undercurrent appears by 
day 30 (not shown) in the opposite direction of the surface 
flow. A cross-section of meridional velocity at 37.4°N on 
day 39 (Figure 8e) shows the core of the poleward 

undercurrent at -500 m depth with a speed of -10 cm/s 
underlying the equatorward surface current which has a 
speed of -20 cm/s. The undercurrent reaches its maximum 
intensity by day 45 (not shown), after which eddies and 
Meddles begin to form and spin off. 

3. Filaments and Eddies 

The upwelling-favorable equatorward winds cause 

cooling of the coastal waters along the coast and produces 
offshore flowing filaments that transport cold water into 
the sea interior. At the surface, several filaments of 

cooler water are generated from upwelling. By day 30, there 
is evidence of filament activity off Cabo da Roca, Cabo de 
Sao Vicente and Cape Ghir (Figure 8a) . The coastal 

irregularities (i.e., capes) off these three coastal 
regions cause surface eddies to spin off from the filaments 
by day 45 (Figure 8b) . By day 60, the formation of some 
smaller mesoscale features in the Gulf of Cadiz is also 
noticeable (Figure 8c) . By day 75, surface eddies appear 
along the entire coastal domain (Figure 8d), with more 
developed features near coastline irregularities, including 
off Cabo da Roca, Cabo de Sao Vicente, Cape Beddouzza and 
Cape Ghir. 

4. Meddles 

At day 60, the development of several subsurface 
anticyclonic eddies off Cape Bedouzza, Cabo de Sao Vicente 
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and Cabo da Roca can be seen at ~1250 m depth (Figure 8f) . 
A cross section of salinity at 39.3°N off Cabo da Roca 
(Figure 8g) shows a salty core with a salinity maximum of 
~35.9 psu at 1250 m depth, which is consistent with the MO 
signature. It has been suggested (Kase et al., 1989) that 
Meddles are generated by the instability of the equatorward 
coastal jet and poleward undercurrent. The deep origin, 
salty signature, and anticyclonic rotation of the eddy west 
of Cabo de Sao Vicente are consistent with Meddy 
observations in this region (Richardson and Tychensky, 
1998) . The Meddy observed off Cabo da Roca (Figure 8f) is 
consistent with results of Kase et al. , (1989), who 
observed eddies in the MO off northern Iberia. 

B. EXPERIMENT TWO 

In Experiment 2, instead of a flat bottom, the 
Gaussian smoothed bottom topography, described in Chapter 
II, was added with the same model forcing conditions (see 
Table 1). In particular, the model was initialized with the 
horizontally averaged annual climatological temperatures 
and salinities. A realistic coastline was used along with 
forcing by the annual climatological winds. 

1. Surface Currents and Upwelling 

In this experiment, the surface current velocities are 
maintained over the shelf rather than close to the shore as 
in the flat bottom case (e.g., compare Figures 8 and 9) . 
The velocities also have the highest magnitudes over the 
continental shelf break, as indicated by the 200 m contour 
in Figure 9. 

A comparison of the surface temperatures indicates 
cooler upwelled water along the entire coastline in 
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Experiment 2 compared to the flat bottom experiment. 
However, the extent of these upwelled waters in Experiment 
2 is more restricted than in Experiment 1, particularly off 
Cape Ghir. This can be explained by the presence of bottom 
topography which traps the flow and opposes the tendency 
for westward propagation due to the planetary beta effect. 
Another area of stronger upwelling is near Eigueira da Eoz, 
where colder surface temperatures are present farther 
offshore in Experiment 2 than in Experiment 1 (e.g., 
compare Eigures 8 and 9). This can be explained as follows: 
At 40.3°N, the 200 m contour depicts a prominent subsurface 
topographic feature which is not present in the flat bottom 
experiment. 

2. Undercurrent 

A poleward undercurrent north of the Gulf of Cadiz 
develops much later in Experiment 2 than in Experiment 1. 
Instead of forming by day 15 as in Experiment 1, the core 
does not form until ~day 33 (not shown) . It reaches its 
maximum intensity by day 45. The core is also deeper in 
Experiment 2 than in Experiment 1 (e.g., compare Eigures 8e 
and 9e) . In Experiment 2, the core of the undercurrent is 
below ~1000 m depth while in Experiment 1 it is shallower 
below ~500 m depth. This is a likely cause for a decrease 
of mesoscale activity in Experiment 2. The deeper 
undercurrent can result in weaker vertical and horizontal 
shears which subsequently reduces the baroclinic and 
barotropic instabilities responsible for the generation of 
mesoscale features (Batteen, 1997; Roed and Shi, 1999). The 
presence of a linearly sloping bottom at the continental 
slope can also produce barotropic shear waves which have 
been demonstrated to be more stable than those found in 
flat bottom cases (Gawarkiewicz, 1991). 
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3. Filaments and Eddies 

Unlike Experiment 1, which had extensive meanders and 
filaments (e.g.. Figure 8c), the only evidence of mesoscale 
phenomena in Experiment 2 is a relatively strong eddy off 

Cabo da Roca by day 45 (Figure 9b) and another weaker eddy 

off Figueira da Foz (Figure 9d) by day 75. Again, the eddy 

off Figueira da Foz, which was not seen in Experiment 1, is 

associated with the topographic feature off the almost 
straight coastline. The filament off Cape Ghir in 
Experiment 1 (Figure 8c) is not present in Experiment 2 
(Figure 9c) . This is likely due to trapping of the surface 
currents by the topography which opposes the westward 
propagation of the upwelled waters. The presence of these 
features demonstrates that bottom topography plays an 
important role in the development and intensification of 
filaments. 

4. Meddles 

Meddy formation occurs in Experiment 2 off Cabo da 
Roca and off Figueira da Foz (Figure 9f) . However, unlike 
Experiment 1, no Meddy forms in Experiment 2 off Cabo de 
Sao Vicente (e.g., compare Figures 8f and 9f) . Note that 
the current follows the smoothed continental shelf in 
Experiment 2 rather than the coastline. In Experiment 2, 
this Gaussian smoothed bottom topography off the capes 
results in the following: Both vortex stretching and the 
increased radius of curvature (compared to the radius in 
Experiment 1) in the Gaussian smoothed topography impede 
the separation of the boundary current. (Note that boundary 
current separation is necessary for the formation of 
Meddles (Martinho, 2003)). As a result, no Meddy forms off 
Cabo de Sao Vicente in Experiment 2. In contrast, off Cabo 
da Roca, there is no significant change in the radius of 
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curvature from Experiment 1 to Experiment 2 (e.g., compare 

Eigures 8f and 9f) . Therefore, Meddies do occur in this 
location for both the flat bottom (Experiment 1) and the 
topography (Experiment 2) cases. 

C. EXPERIMENT THREE 

Experiment 3 uses the same Gaussian smoothed bottom 
topography as in Experiment 2. However, this experiment 
uses the full annual temperature and salinity climatology, 
which models the observed density/sea height gradient and 
introduces thermohaline forcing into the model. The same 
forcing by annual climatological winds and a realistic 
coastline was used. 

1. Surface Currents and Upwelling 

The equatorward surface current in Experiment 3 is 
comparable to that in Experiment 2 (e.g., compare Eigures 9 

and 10) . The coastal current is slightly weaker in 
Experiment 3 due to the thermohaline forcing which opposes 
the wind forcing. Geostrophic adjustment can produce a 
velocity smaller in magnitude but in opposition to the main 
flow which reduces the current speed. 

The difference in the coastal area temperatures in 
Experiments 2 and 3 reflect the different climatologies 
used in the model. Eor example, to the north of Eigueira da 
Eoz, there are differences in both the temperatures and the 
upwelling extent (e.g., compare Eirgures 9a and 10b). 

2. Undercurrent 

A comparison of cross-sections of the meridional 
velocities for Experiments 2 and 3 (Eigures 9e and lOf) 
shows a shallower extent of the undercurrent in Experiment 
3, but with a well defined core of similar magnitude at 
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~1000 m depth. As in Experiment 2, the deep location of the 
undercurrent results in weaker vertical and horizontal 
shears. This reduces the baroclinic and barotropic 
instabilities that generate mesoscale features. 

3. Filaments and Eddies 

In Experiment 3, broad filaments (e.g., Eigure 10b) 
occur to the north off Cabo da Roca. A meander present off 
Cabo da Roca in Experiment 3 (e.g., Eigure 10c) is 

comparable to that in Experiment 2, except that the 
intensity of the surface eddy offshore of the filament is 
much less than in Experiment 2 (e.g., compare Eigures 9c 

and lOd). Due to the effect of full climatology used, which 
tightens the surface current near the coast, the extent of 
the filament off Cape Ghir is also less in Experiment 3 
(e.g., compare Eigures 9d and lOe) . In Experiment 3, no 
surface eddy spins off Eigueira da Eoz as was observed in 
Experiment 2 (e.g., compare Eigures 9d and lOe). 

4. Iberian Current 

In Experiment 3, thermohaline forcing was the only 
addition incorporated. A new feature in Experiment 3 is the 
formation of a poleward flow north of ~39°N (centered at 
~10°W) offshore of the coastal equatorward flow which is 
discernible by day 5 (Eigure 10a) . This is consistent with 
the IC, which flows off the IP. 

Other models have produced a similar poleward flow 
such as Eigure 15 of Coelho et al. (2002), which is shown 
as Eigure 11. A comparison of Eigures 10a and 11 reveals 
that both model runs have an equatorward surface current 
nearshore and an offshore poleward current between ~39°N 
and ~42°N centered at ~10°W. Coelho et al. (2002) suggested 

a number of explanations for the surface poleward flow 
including: the type of wind forcing used (i.e.. 
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climatological versus actual (in their case, 1994) winds); 
remnant of a winter surface poleward current (generated by 
thermohaline forcing); and wind stress curl. In Experiment 
2, no surface poleward flow was generated in this region. 
As a result, we conclude that of all the proposed 
mechanisms by Coelho et al., (2002), the surface poleward 
flow is likely generated by the additional thermohaline 
forcing instead of either by wind forcing or wind stress 
curl. 

5. Meddles 

Meddy formation occurs in Experiment 3 first off Cabo 
da Roca and then off Eigueira da Eoz (Eigures lOg and lOh). 
As in Experiment 2, no Meddy forms off Cabo de Sao Vicente. 

A comparison of the results of all three experiments 
shows that while most features are simulated in all the 
experiments, unexpectedly a realistic subsurface mesoscale 
feature (i.e., a Meddy off Cabo de Sao Vicente) only occurs 
in the flat bottom NCCS model (Experiment 1) . In the next 
section, we will investigate if this is a consequence of 
the type of smoothed bottom topography used in both 
Experiments 2 and 3. 

D. EXPERIMENT FOUR 

Instead of Gaussian smoothed bottom topography, direct 
iterative topography, described earlier, is used in this 
Experiment. Note that in the direct iterative smoothing 
method, the continental shelf more closely approximates the 
raw topography than the Gaussian smoothing method (e.g., 
compare Eigures Ib-ld and also the 200 m isobath in Eigures 
10 and 12) . The same model forcing conditions (see Table 1) 
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are used as in Experiment 3. Specif ically, the model is 
initialized with the full annual climatological 
temperatures and salinities and annual wind forcing. 

1. Surface Currents and Upwelling 

By day 5, a weak IC develops north of 39°N but is not 
subsequently maintained (likely due to the incorporation of 
annual rather than winter thermohaline forcing). By day 30, 
an equatorward surface current develops along the 
continental shelf break and upper continental slope (Figure 
12b) . Strong upwelling also occurs off the west coast of 
Morocco due to strong equatorward winds in this region. By 
day 45 (Figure 12c) , upwelling is enhanced near the capes 
off the IP. By day 60, upwelling generally has a greater 
westward extent in Experiment 4 than in Experiment 3 (e.g., 

compare Figures lOd and 12d) . By day 75 (Figure 12e) , the 
offshore extent of upwelled water off Cabo da Roca and 
Figueira da Foz extend well seaward of the continental 
shelf break. 

2. Undercurrent 

In Figure 12f, the results for day 39 of the 
meridional velocity are shown in a cross-section for 
Experiment 4. Here, the poleward undercurrent is more 
confined to the slope than in Experiment 3 (e.g., compare 

Figures lOf and 12f). The core of the poleward undercurrent 
is located at ~800 m depth in Experiment 4 rather than at 
~1000 m depth in Experiment 3. This shallower core of the 
undercurrent increases the vertical and horizontal shears 
and is likely responsible for the generation of greater 
mesoscale activity off Cabo de Sao Vicente (e.g., compare 
Figures lOe and 12e). 
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3. Filaments and Surface Eddies 

In the previous Gaussian topography experiment, a weak 
filament developed off Cape Ghir by day 75 (Figure 12e). In 
Experiment 4, a filament develops off Cape Ghir much 
earlier, i.e., by day 30 (Figure 12b), and subsequently 
spins off a surface eddy (day 45, Figure 12c) . A filament 
not found in Experiment 3 off Cabo de Sao Vicente is now 
present in Experiment 4 (e.g., compare Figures lOe and 

12e) . In particular, a surface eddy develops to the south 
of Cabo de Sao Vicente in the Gulf of Cadiz. While 
filaments extend offshore from both Cabo da Roca and 
Figueira da Foz at the same extent as Experiment 3 by day 
45 (e.g., compare Figures 10c and 12c), a surface eddy does 

not spin off Cabo da Roca (i.e., by day 75, Figure 12e, 
compared to day 60, Figure lOd) . Another feature not found 
in Experiment 3 but now found in Experiment 4 is surface 
eddy activity in the Gulf of Cadiz (e.g., compare Figures 
lOe and 12e) . 

4. Meddles 

A Meddy develops to the west of Cabo de Sao Vicente by 
day 30 (Figure 12g), which is consistent with available 
observations of Meddles (e.g., Batteen et al., 2000). This 

Meddy was also produced in the flat bottom case (e.g.. 
Figure 8f) , but not in both Gaussian smoothed topography 
cases. Experiments 2 and 3 (e.g.. Figures 9f and lOg) . Just 

as in Experiment 2, a Meddy forms off Figueira da Foz (not 
shown), although much weaker in intensity (compare Figures 
lOg and 12e). As in previous experiments, a Meddy off Cabo 
da Roca forms (e.g.. Figure 12e) and subsequently 
propagates to the west (not shown) . 
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V. EFFECTS OF NORTHERN CANARY CURRENT SYSTEM 
FEATURES ON SOUND SPEED STRUCTURE 

A. INTRODUCTION 

The model experiments in the previous chapter 
reproduced the major characteristics and unique features of 
the NCCS and highlighted the observed mesoscale variability 
in the NCCS. This section demonstrates that these features 
lead to considerable variability in sound speed which can 
have a significant impact on sound speed structure. 

The sound speed structure in a region is affected by 
two general processes: vertical displacement of ocean 
layers and horizontal advection of water masses (Lurton, 
2002) . Ocean layers are displaced vertically by features 
such as internal waves, and surface and subsurface eddies, 
which create variations in water density with depth. 
Upwelling, the lifting of deep water masses to the upper 

ocean, also vertically displaces water. Horizontal 
advection occurs when meanders and filaments in the surface 
current carry warmer water seaward, transporting large 

volumes of warm water. 

Applying this to the NCCS region, we note that the 
NCCS is affected by upwelling along the coast which 
displaces water vertically. On the horizontal scale, the 
Canary Current advects cool fresh water equatorward while 
the undercurrent advects warm saline water poleward. There 
also is the advection of MO by the undercurrent which is 
unique to the NCCS region. Filaments can carry cooler water 
offshore, while eddies that detach from the current near 
the capes can transport large volumes of cool water. In 

addition, all of these processes that change the sound 
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speed are variable in time (as shown in Chapter IV) which 
can further complicate the acoustic picture. 

In this chapter, we examine the sound speed structure. 
In particular, sound speed is calculated and plotted for 
horizontal levels, such as the surface. Sound speed and 
sound speed change are also calculated and plotted in 
vertical cross sections. These are used to determine how 
MO, upwelling, undercurrents, filaments, surface eddies and 
Meddles change the sound speed. Sound speed is calculated 
at the time of analysis and the change of sound speed is 
the difference in sound speed from the initial time 
(climatology) to the analysis time. 

B. FEATURES 

1. Mediterranean Outflow (MO) 

MO is generally found throughout the NCOS model 
domain. For example, a cross-section of salinity through 
the Gulf of Cadiz at 36°N (Figure 13a) shows the MO 
signature centered at ~1200 m depth and also the 
characteristically high salinity values of the North 
Atlantic waters in the upper 300 m depth. To the south 
(32°N, Figure 13b) and to the north (39°N, Figure 13c), 
evidence of MO centered at ~1200 m depth can also be seen. 

For the same day and latitudes, the MO signature is 
discernible in the sound speed profiles (Figures 13d, 13e, 

and 13f) creating two distinct sound channels centered at 
~500 m and ~2000 m depths. This particular profile is 
unique to the NCOS region due to warm, salty waters of the 
Mediterranean mixing with the North Atlantic. Important 
acoustically, the double minimum features are strong not 
only in the Gulf of Cadiz (Figure 13d) but also towards the 
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south (Figure 13e) and north (Figure 13f), exemplifying the 
same spreading of the MO seen in salinity cross sections 
(Figures 13a, 13b and 13c) . 

2. Upwelling 

In Experiment 4, a band of upwelled water is visible 
along the shore at the surface (Figure 12e) . As expected, 
the surface sound speed closely resembles the temperature 
plot with the slowest sound speeds in areas where the 
temperatures are lower (e.g., compare Figures 12e and 14a) . 
For example, a cross-section at 33.5°N (Figure 14b) shows 
that the upward sloping isotherms over the shelf are 
consistent with the presence of upwelling. There is also an 
upward slope in sound speed with the slower speeds in the 
areas of cooler waters (Figure 14c). The significantly 
cooler upwelled water subsequently results in a decrease of 
~25 m/s in sound speed from the surface to ~500 m depth 
near the shelf (Figure 14d) . 

In the region of upwelling, sound speed fluctuation is 
driven by the change in temperature. Although upwelling 
primarily impacts the mixed layer, the horizontal gradient 
of sound speed in the upper layer is not widespread and 
changes ~25 m/s over a ~150 km distance. 

3. Filaments and Surface Eddies 

Several eddies of cooler water are generated on the 
surface. In particular, by day 75 in Experiment 4 (Eigure 
12e) , a well-defined, cold-core eddy centered at ~11°W has 
spun off from these filaments at Cape Ghir. The surface 
sound speed is lower in the eddy than in the surrounding 
waters (Eigure 14a) . The lower sound speed in the eddy can 
extend to over 200 m depth. The decrease in sound speed 
from initial time is ~20 m/s across a ~200 km horizontal 
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distance (Figure 14d) . Note that the sound speed decrease 
is lowest in the core of the eddy. 

These mesoscale surface eddies are significant 
acoustical features due to their large gradients of sound 
speed. The sound speed structure in the eddy changes over a 
relatively short longitudinal distance, leading to 
relatively large gradients. 

4. Meddles 

Meddles developed in all the experiments off Cabo de 
Roca. A typical cross-section of salinity shows that the 
core of the Meddy (e.g.. Figure 15a) is salty which 
corresponds to a higher sound speed (Figure 15b) . At 39.7°N 
off Cabo da Roca, a cross-section of sound speed change 
(Figure 15c) shows that the salty core at 1250 m depth 
(Figure 15d) increases the sound speed. In particular, the 
sound speed across the Meddy increases by ~6 m/s over a 
~100 km horizontal and ~1000 m vertical region (Figure 
15c) . Whereas temperature is a foremost factor affecting 
sound speed, the NCCS represents a region uniquely affected 
by the highly saline MO. 

5. Undercurrent 

The undercurrent (e.g.. Figure lOf) advects 

subtropical water poleward along the western coasts. North 
of the Gulf of Cadiz, the undercurrent advects MO waters 
and has a significant effect on sound speed structure. For 
example, at 37.4°N (Figure 16a), the downward sloping 
isotherms near the shelf at ~1000 m depth and below are 
consistent with downwelling (Figure 16a) in conjunction 
with the subsurface advection of subtropical water. A cross 
section of sound speed difference at 37.4°N (Figure 16b) 
shows a ~6 m/s increase due the added salinity of the MO to 


34 



the warm undercurrent. This increase in sound speed occurs 
over a ~50 km horizontal and a ~500 m vertical extent. 
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VI. SUMMARY 


Several process-oriented experiments were conducted to 
determine if different types of topographic smoothing, 
e.g., Gaussian and direct iterative methods, could 
significantly influence the generation, evolution, and 
maintenance of the classical as well as unique features of 
the NCOS. The incorporation of a volume constraint to the 
terrain-following numerical ocean model, in this case the 
POM, was shown to improve model stability and the overall 
performance on the long term. 

While most of the NCOS features were realistically 
simulated, a key comparison of the results showed that 
unexpectedly a realistic subsurface mesoscale feature (i.e. 
a Meddy off Cabo de Sao Vicente) was simulated in both the 
flat bottom and iterative bottom topography NCOS models but 
not in the model using Gaussian bottom topography. The 
iterative bottom topography was shown to better preserve 
the raw topography and subsequently produced the subsurface 
feature. The results of these process-oriented experiments 
emphasize that numerical models of ocean circulation 
require important choices, which are both numerical and 
physical. 

The effects of these mesoscale features in the current 
system on the acoustic structure in the region were 
systematically examined. Recall in section V.B.l that to 
the north of the Gulf of Cadiz, the added salinity of the 
MO caused a significant increase in sound speed (e.g., 
compare Figure 13 a-c to Figures 13d-f) . 

Table 2 summarizes the other effects of sound speed of 
each feature discussed in the previous section (i.e.. 
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Chapter V) for this region. Results are presented in terms 
of effective change in sound speed from initial conditions 
as well as spatial extent. Of critical acoustic importance 
are the spatial gradients of sound speed caused by these 
features. 


NCCS 

Feature 

Ac 

Lx 

Lz 

Upwelling 

-25 m/s 

150 km 

500 m 

Surface Eddies 

-20 m/s 

2 0 0 km 

2 0 0 m 

Meddies 

+ 6 m/s 

5 0 km 

1000 m 

Undercurrent 

+ 6 m/s 

5 0 km 

500 m 


Table 2. Impact of features on sound in the 

NCOS . 

In the NCCS, upwelling clearly experiences the 
greatest change on sound speed, affecting a large surface 
area through the mixed layer. Surface eddies which spin off 
upwelling filaments also cause a relatively large decrease 
in sound speed. The Meddles and undercurrent presented the 
most intense horizontal and vertical gradients of sound 
speed change. 

Overall, this study demonstrates that features of 
classic coastal/littoral eastern boundary current systems, 
like the NCCS, do have a significant impact on the sound 
speed structure. Knowledge gained from this study should 
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allow Navy Oceanographers to better understand the complex 
ocean dynamics observed in the littoral/coastal systems and 
more accurately predict sound speed variability in these 
regions. In turn, the Navy will be better equipped to 
support mine and anti-submarine warfare operations in these 
areas. 
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APPENDIX 


This appendix contains all 


figures referred to 


throughout the text. 
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-12 -10 -8 -6 
Longitude 

Figure la. Model domain for the Northern Canary Current 
System (NCCS) is bounded by 31.5°N to 41.5°N, 16.5°W (not 

shown) to 6°W. Geographical locations are labeled. 


41 



Latitude 


Raw Topography 



1 3000 

2000 

1000 

-0 

-1000 

-2000 

-3000 

--4000 



-16 -14 -12 -10 -8 

Longitude 


-6 


Figure lb. Raw topography, 
Sandwell and Smith (1996) . 
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Figure Ic. Smoothed topography, depths in m, obtained from 
Sandwell and Smith (1996) after applying a linear two- 
dimensional low-pass filter. 
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Figure Id. Iterative topography, depths in m, obtained 
Sandwell and Smith (1996) after applying a linear 
dimensional low-pass filter. 
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Annual Surface Temperature 
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Figure 2. Levitus annual climatological surface temperature 
(°C) . 
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Annual Surface Salinity 
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Figure 3. Levitus annual climatological surface salinity 
(°C) . 
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Annual Wind Forcing 



Figure 4. NCCS wind stress vectors (Pascals) calculated 
from annual climatological ECMWF winds (From: Trenberth et 
al., 1990) 
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Figure 5a. Raw topography for the Northern Canary Current 
System, depths in meters. Contour lines at 100, 200, 500 
and 1000 m depth. 
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Figure 5b. Topography smoothed 
dimensional filter method for the 
System, depths in meters. Contour 
and 1000 m depth. 
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Figure 5c. Topography smoothed with the direct iterative 
method for the Northern Canary Current System, depths in 
meters. Contour lines at 100, 200, 500 and 1000 m depth. 
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Figure 5d. Cross sections of topography across a seamount 
at 36.4°N showing the methods of topography smoothing. The 
blue line is the raw (unsmoothed topography). The red line 
is the topography smoothed with a Gaussian filter. The 
black line is the topography smoothed with the alternative 
POM method. The dashed magenta line is the topography 
smoothed using the one-dimensional direct iterative 
technique. 
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w/ Volume Constraint, Day 1, Surface. 
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Figure 6a. Surface temperature (°C) and velocity vectors 
(cm/s) for day 1 with Volume Constraint. 
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w/ Volume Constraint, Day 90, Surface. 
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Figure 6b. Surface temperature (°C) and velocity vectors 
(cm/s) for day 90 with Volume Constraint. 
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w/o Volume Constraint, Day 1, Surface. 
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Figure 7a. Surface temperature (°C) and velocity vectors 
(cm/s) for day 1 without Volume Constraint. 
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w/o Volume Constraint, Day 90, Surface. 
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Figure 7b. Surface temperature (°C) and velocity vectors 
(cm/s) for day 90 without Volume Constraint. 
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EXP 1, Day 30, Surface. 
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Figure 8a. Surface temperature (°C) and velocity vectors 
(cm/s) on day 30 for Experiment 1. 
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EXP 1, Day 45, Surface. 
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Figure 8b. Surface temperature (°C) 
(cm/s) on day 45 for Experiment 1. 
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EXP 1, Day 60, Surface. 
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Figure 8c. Surface temperature (°C) and velocity vectors 
(cm/s) on day 60 for Experiment 1. 
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EXP 1, Day 75, Surface. 
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Figure 8d. Surface temperature (°C) 
(cm/s) on day 60 for Experiment 1. 
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EXP 1, Velocity Cross-section at 37.4N, Day39. 




Figure 8e. Velocity (m/s) cross-section at 37.4°N on day 
for Experiment 1. 
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EXP 1, Day 60,1250m Depth. 
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Figure 8f. Salinity (psu) and velocity vectors (cm/s) 
1250 m depth on day 60 for Experiment 1. 
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EXP 1, Salinity Cross-section at 39.3N, Day60. 
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Figure 8g. Salinity (psu) cross-section at 39.3°N on day 
for Experiment 1. 
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EXP 2, Day 30, Surface. 
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Figure 9a. Surface temperature (°C) and velocity vectors 
(cm/s) on day 30 for Experiment 2. 
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EXP 2, Day 45, Surface. 
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Figure 9b. Surface temperature (°C) 
(cm/s) on day 45 for Experiment 2. 
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EXP 2, Day 60, Surface. 
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Figure 9c. Surface temperature (°C) 
(cm/s) on day 60 for Experiment 2. 
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EXP 2, Day 75, Surface. 
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Figure 
(cm/s) 


9d. Surface temperature 
on day 75 for Experiment 2 
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EXP 2, Velocity Cross-section at 37.4N, Day39. 
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Figure 9e. Velocity (m/s) cross-section at 37.4°N on day 
for Experiment 2. 
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EXP 2, Day 60,1250m Depth. 
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Figure 9f. Salinity (psu) and velocity vectors (cm/s) 
1250 m depth on day 60 for Experiment 2. 
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EXP 3 , Day 5, Surface. 
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Figure 10a. Surface temperature (°C) and velocity vectors 
(cm/s) on day 5 for Experiment 3. 
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EXP 3 , Day 30, Surface. 
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Figure 10b. Surface temperature (°C) and velocity vectors 
(cm/s) on day 30 for Experiment 3. 
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EXP 3 , Day 45, Surface 
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Figure 10c. Surface temperature (°C) and velocity vectors 
(cm/s) on day 45 for Experiment 3. 
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lOd. Surface temperature (°C) and velocity vectors 
on day 60 for Experiment 3. 
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EXP 3 , Day 75, Surface. 
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Figure lOe. Surface temperature (°C) and velocity vectors 
(cm/s) on day 75 for Experiment 3. 
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EXP 3, Velocity Cross-section at 37.4N, Day39. 




Figure lOf. Velocity (m/s) cross-section at 37.4°N on 
39 for Experiment 3. 
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EXP 3, Day 75,1250m Depth. 




Figure lOg. Salinity (psu) and velocity vectors (cm/s) 
1250 m depth on day 75 for Experiment 3. 
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Figure 11. Zoom of monthly averaged velocity field at level 
1 (5 m) for July 1994. To avoid clutter, the velocity 

vectors are plotted every second grid point in both 
directions (from Coelho et al. (2002) Figure 15). 
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EXP 4, Day 5, Surface. 
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Figure 12a. Surface temperature (°C) 
(cm/s) on day 5 for Experiment 4. 
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EXP 4, Day 30, Surface. 
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Figure 12b. Surface temperature (°C) and velocity vectors 
(cm/s) on day 30 for Experiment 4. 
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12c. Surface temperature (°C) and velocity vectors 
on day 45 for Experiment 4. 
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EXP 4, Day 60, Surface. 
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Figure 12d. Surface temperature (°C) and velocity vectors 
(cm/s) on day 60 for Experiment 4. 
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Figure 12e. Surface temperature (°C) and velocity vectors 
(cm/s) on day 75 for Experiment 4. 
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EXP 4, Velocity Cross-section at 37.4N, Day39. 




Figure 12f. Velocity (m/s) cross-section at 37.4°N on 
39 for Experiment 4. 
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EXP 4, Day 30,1250m Depth. 
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Figure 12g. Salinity (psu) and velocity vectors (cm/s) 
1250 m depth on day 30 for Experiment 4. 
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EXP 4, Salinity Cross-section at 36N, DayO. 
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Figure 13a. Salinity (psu) cross-section at 36°N on day 0 
for Experiment 4. 
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EXP 4, Salinity Cross-section at 32N, DayO. 
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Figure 13b. Salinity (psu) cross-section at 32°N on day 0 
for Experiment 4. 
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EXP 4, Salinity Cross-section at 39N, DayO. 



Figure 13c. Salinity (psu) cross-section at 39°N on day 
for Experiment 4. 
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EXP 4, C Cross-section at 36N, Day 0. 



Figure 13d. Sound speed (m/s) cross-section at 36°N on 
0 for Experiment 4. 
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EXP 4, C Cross-section at 32N, Day 0. 
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Figure 13e. Sound speed (m/s) cross-section at 32°N on day 
0 for Experiment 4. 


































EXP 4, C Cross-section at 39N, Day 0. 
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Figure 13f. Sound speed (m/s) cross-section at 39°N on day 
0 for Experiment 4. 
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EXP 4, Day 75, Surface. 
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Figure 14a. Surface sound speed (m/s) and velocity vectors 
(cm/s) on day 75 for Experiment 4. 
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EXP 4, Temp Cross-section at 33.5N, Day75. 
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Figure 14b. Temperature (°C) cross-section at 33.5°N on day 
75 for Experiment 4. 
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EXP 4, C Cross-section at 33.5N, Day 75. 



Figure 14c. Sound speed (m/s) cross-section at 33.5°N 
day 75 for Experiment 4. 
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EXP 4, Delta C Cross-section at 33.5N, Day75. 



Figure 14d. Difference in sound speed 
conditions, cross-section at 33.5°N 
Experiment 4. 
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EXP 3 , Salinity Cross-section at 39.7N, Day75. 
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Figure 15a. Salinity (psu) cross-section at 39.7°N on 
75 for Experiment 3. 
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EXP 3, C Cross-section at 39.7N, Day 75. 
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Figure 15b. Sound speed (m/s) cross-section at 39.7°N 
day 75 for Experiment 3. 



on 


95 












EXP 3, Delta C Cross-section at 39.7N, Day75. 




Figure 15c. Difference in sound speed 
conditions, cross-section at 39.7°N 
Experiment 3. 
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EXP 3, Day 75,1250m Depth. 
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15d. Salinity (psu) and velocity vectors (cm/s) 
depth on day 75 for Experiment 3. 
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Figure 16a. Temperature (°C) cross-section at 37.4°N on day 
39 for Experiment 3. 
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EXP 3, Delta C Cross-section at 37.4N, Day39. 




Figure 16b. Difference in sound speed 
conditions, cross-section at 37.4°N 
Experiment 3. 
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